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Abstract: The BL Lac object (blazar) Mrk 421 was observed during its outburst in April 2004 with the 
Whipple 10 m telescope for a total of about 24.5 hours. The measured 7-ray rate varied substantially 
over the range from 4 to 10 7's/min and eventually exceeded the steady 7-ray rate of the Crab Nebula 
(standard candle) by a factor of 3. The overall significance of the 7-ray signal exceeded 70 a and the 
total number of excess events was more than 10,000. The signal light curve does not show any particular 
variability pattern. This unique Mrk 421 outburst enabled the measurement of a high quality spectrum of 
very high-energy 7 rays in a high state of emission. This spectrum is a power-law and it extends beyond 
10 TeV. 
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Introduction 

Mkn 421 is the first detected, the closest known 
(red shift z = 0.030), and one of the best studied 
TeV 7-ray emitting blazars. The very high-energy 
(VHE) 7-ray emission arises from the particles ac- 
celerated in a relativistic jet directed along our line- 
of-sight. Since its discovery, Mkn 421 has shown 
very low baseline TeV 7-ray emission with a few 
extremely rapid flares on time-scales from one day 
to 15 minutes [7]. In 2000 and 2001 Mkn 421 went 
into a flare state with an average flux of 4 times that 
of the Crab Nebula. Data taken during this flare 
have been used to extract the energy spectrum at 
high energies, up to 20 TeV. Noticeable variations 
in the hardness of the TeV 7-ray energy spectrum 
have been reported [1, 12]. A number of success- 
ful multi-wavelength campaigns for Mkn 421 have 
revealed evidence of a correlation of the simulta- 
neously measured fluxes in X-ray and TeV 7-ray 
energy bands [11,4]. 

A one-zone synchrotron-self-Compton (SSC) 
model described the observational results rea- 
sonably well (e.g., [10]). During April 2004, an 
intensive multi-wavelength monitoring campaign 
on the TeV blazar Mkn 421 has been performed 
simultaneously in radio, optical, X-ray, and 7 
rays. The source was seen to be active in X-rays 
and TeV 7 rays with the peak flux exceeding 



the 3 Crab level [4]. The time-averaged energy 
spectrum of Mkn 421 during the flaring state has 
been measured at high energies with H.E.S.S. 
using large zenith angles observations [2], and 
MAGIC [3]. The Whipple 10 m telescope was 
observing Mkn 421 extensively during the April 
2004 flare. Here we report the results on the TeV 
7-ray energy spectrum of Mkn 42 1 measured in a 
high flaring state during the April 2004 outburst 
with the Whipple 10 m telescope. 

The 10 m Whipple Telescope 

The reflector of the Whipple Observatory imaging 
atmospheric Cherenkov telescope (IACT) (see Fig- 
ure 1) is a tessellated structure consisting of 248 
spherical mirrors, which are hexagonal in shape 
and 61 cm from apex to apex, arranged in a hexag- 
onal pattern. The mirrors are mounted on a steel 
support structure, which has a 7.3 m radius of cur- 
vature with a 10 m aperture. Each individual mir- 
ror has — 14.6 m radius of curvature and is pointed 
toward a position along the optical axis at 14.6 m 
from the reflector. This arrangement constitutes a 
Davies, Cotton design [5] of the optical reflector. 
The point-spread function of the telescope has a 
FWHM of - 7.2' on-axis. 

In 1999, a high-resolution camera (GRANITE III) 
was installed at the telescope [6]. It consists 
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RA: 110426.3 DEC: 381237.7 J2000.0 



Figure 1: Whipple Observatory 10 m Telescope at 
Mount Hopkins, Arizona. 

of 379 photo-multiplier tubes (PMTs) in a close 
packed hexagonal arrangement (each PMT sub- 
tending 0.11° on the sky) and has a 2.6° diameter. 
A set of light concentrators is mounted in front of 
the pixels to increase the light-collection efficiency 
by —38%. The camera triggers if the signal in each 
of at least 3 neighboring PMTs out of the inner 
331 exceeds a threshold of 32 mV, corresponding 
to -8-10 photo-electrons. The post-GRANITE III 
upgrade trigger rate of the Whipple Observatory 
10 m telescope is —20-30 Hz at zenith. 

Observations & Analysis 

Mkn 421 was observed with the Whipple Observa- 
tory 10 m IACT for about 24.5 hours of on-source 
data between April 9th and April 23d, 2004. Data 
have been taken in good weather at zenith angles 
less than 30° (mean elevation is about 75°). Dur- 
ing observations the raw detection rate was about 
22 Hz. Observations have been performed using 
tracking mode in data runs of 28 minutes each. The 
recorded images are first flat-fielded using nightly 
measured nitrogen arc lamp PMT responses and 
then cleaned by applying a standard picture and 
boundary technique with canonical thresholds of 
4.25 and 2.25 times the standard deviation of the 
PMT pedestal distributions, respectively (see, e.g., 
[9]). To characterize the shape and orientation 
of calibrated images, the standard second-moment 
image parameters are calculated as described in 
[16]. The CAnalyze package developed at Pur- 
due University [13] is used for the primary data 
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Figure 2: The two-dimensional map of excess sig- 
nificances over the sky region around Mkn 421 
generated for the data sample corresponding to 
a high emission state. The map of uncorrected 
rectangular angular bins of a 0.1° size has been 
smoothed using Gaussian kernel with the angular 
size of 0.22°. The maximal significance of 7-ray 
excess using this analysis is 12.3 a. 

Table 1 : Supercuts selection criteria for the Whip- 
ple Observatory 10 m telescope (dc stands for dig- 
ital counts). 



Quantity 


Image Parameter Cut 


Trigger 


1st & 2nd brightest pixel > 30 dc 


Shape 


0.05° < width < 0.12° 




0.13° < length < 0.25° 


Muon cut 


length/ size < 0.0004 (° dc" 1 ) 


Quality cut 


0.4° < distance < 1.0° 


Orientation 


a < 15° 



analysis. Analysis methods known as Supercuts 
(see Table 1) and extended Supercuts are applied 
for the 7-ray/hadron separation. The latter uti- 
lizes dynamic, energy dependent orientation and 
shape cuts. The significance map of the sky re- 
gion around Mkn 421 for the data set correspond- 
ing to a high emission state (exposure of 55 min- 
utes) is shown in Figure 3. The energy spectrum 
of Mkn 421 is reconstructed using algorithms de- 
scribed in [15]. 

Simulations 

The KASCADE shower simulation code [8] is 
used for generating the 7-ray and cosmic -ray in- 
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Figure 3: Light curve of TeV 7-ray emission from 
Mkn 421 during its outburst in April 2004. 
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duced air showers within the corresponding range 
of zenith angles and in the primary energy range 
between 50 GeV and 100 TeV, assuming the 7-ray 
energy spectrum to be a power-law of 2.5. Sim- 
ulations of the response of the Whipple Observa- 
tory 10 m telescope are carried out using the GrISU 
code, developed by the Grinnell College and Iowa 
State University groups. 

Results 

The light curve of Mkn 421, measured with the 
Whipple Observatory 10 m telescope during the 
outburst in April 2004, does not reveal any particu- 
lar pattern (see Figure 2). Two data runs with a flux 
above 3 Crab and an exposure of 55 minutes give 
exceptionally high 7-ray rate with an average value 
of 9.85 min -1 and a total of 545 excess counts. 
This excess corresponds to a significance of the 7- 
ray signal at a level of 16.4 a. The correspond- 
ing 7-ray energy spectrum is consistent with a pure 
power-law of index 2.66. This spectrum does not 
indicate any apparent gradual change of slope (see 
Figure 4). Such unusual behavior of the Mkn 421 
TeV 7-ray spectrum measured in a high emission 
state suggests that the previously observed curva- 
ture in Mkn 421 7-ray spectrum is an intrinsic fea- 
ture of the source rather than the footprint of 7-ray 
absorption by interacting with the soft photons of 
extragalactic background light. 

Modeling 

To model the multi-wavelength spectrum of the 
BL Lac object Mkn 421 in a homogeneous Syn- 



Figure 4: The energy spectrum of Mkn 421 mea- 
sured in a high emission state during the outburst 
in April 2004. 

chrotron Self Compton (SSC) scenario we use an 
approach described in [14]. This method involves 
prescribing an injection function for relativistic 
electrons and solving the two time-dependent ki- 
netic equations for the electron and photon distri- 
butions of the source. All relevant physical pro- 
cesses are taken into account in the code, i.e., 
synchrotron radiation, inverse Compton scattering, 
photon-photon pair production, and synchrotron 
self-absorption. Seven model parameters are re- 
quired to specify a source in a stationary state. 
These are the Doppler factor 6, the radius R of the 
source, the magnetic field strength B, the index s 
of the electron injection spectrum, the Lorentz fac- 
tor at the cut-off of the injection spectrum j m ax, 
the amplitude of the injection spectrum q e , and the 
effective escape time of relativistic electrons t esc . 
To optimize a fit to a particular data set, we used 
physically motivated starting values for these seven 
parameters (see [10]). This is done by identifying 
six scalars characterizing the blazar spectrum, the 
peak frequency v sync of the synchrotron emission, 
the peak frequency v com pt of the inverse Compton 
emission, the total nonthermal luminosity L, the 
approximate ratio r\ of the total flux in the inverse 
Compton part to that in the synchrotron part of 
the spectrum, the break frequency inbreak, the low- 
frequency spectral index, and finally the fastest 
variability timescale t var . These observables en- 
able one to find reasonable starting values of the 
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Table 2: Summary of Mkn 421 data analyzed by standard supercuts criteria (see Table 1.) 





Time [min] 


ON 


OFF 


Excess 


i? 7 min 1 


S/N [a] 


April 2004 


1472 


18383 


6872 


11511 


7.82 


71.82 


High state 


55 


825 


280 


545 


9.85 


16.40 
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Figure 5: The spectral energy distribution of the 
Mkn 421 emission. The X-ray data, correspond- 
ing to low, medium, and high emission states have 
been taken from [4] . 

seven parameters of the SSC model. Using this ap- 
proach we tried to fit the measured spectral energy 
distribution for a high state of Mkn 42 1 . Despite 
that the X-ray spectra can be reasonably well re- 
produced by the one-zone SSC model, it does not 
provide a satisfactory fit to the TeV 7-ray spectrum 
of Mkn 421 measured in a high emission state (see 
Figure 5). It is worth noting that unfolding the TeV 
7-ray spectrum from the IR absorption (see [ 1 0] for 
further details) yields an energy spectrum substan- 
tially harder that the measured one, which is more 
difficult to interpret in the framework of the one- 
zone SSC model. 

Discussion 

Observations of Mkn 421 with Whipple [12] and 
HEGRA [1] during its historical flare in 2000 and 
2001 have been used to extract the energy spec- 
trum at high energies, up to 20 TeV. As stated by 
both groups, the energy spectrum of Mkn 421 is 
evidently curved. Analysis of the data revealed 
significant variations of the spectral slope at ener- 
gies below 3 TeV depending on the emission state. 
However the best empirical fit to the data taken in 
both the high and the low states gives the same cut- 



off energy of about 3.6 TeV. Recently, the gradual 
steepening in the Mkn 421 TeV energy spectrum 
was evidently detected with HESS [2] and Magic 
[3] during their 2004 and 2005 observational cam- 
paigns. This is consistent with the assumption that 
all Mkn 421 spectra are affected by intergalactic 
absorption. Here we report on the TeV 7-ray spec- 
trum as measured with the Whipple 10 m telescope 
during a high emission state (at the flux level of 3 
Crab), which does not indicate a cut-off feature in 
muti-TeV energy range and can be well-fitted by 
a pure power-law. Such an unusual spectrum in- 
dicates that the curvature could be an intrinsic fea- 
ture of the source. This spectrum offers a challenge 
for the one-zone SSC model, and it may severely 
constrain the spectral energy distribution of extra- 
galactic background light, which is the subject of a 
forthcoming refereed publication. 
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